Nonlinear heat transport in mesoscopic conductors: Rectification, Peltier effect and 

Wiedemann-Franz law 
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We investigate nonlinear heat properties in mesoscopic conductors using a scattering theory of transport. 
Our approach is based on a leading-order expansion in both the electrical and thermal driving forces. Beyond 
linear response, the transport coefficients are functions of the nonequilibrium screening potential that builds up 
in the system due to interactions. Within a mean-field approximation, we self-consistently calculate the heat 
rectification properties of a quantum dot attached to two terminals. We discuss nonlinear contributions to the 
Peltier effect and find departures from the Wiedemann-Franz law in the nonlinear regime of transport. 
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I. INTRODUCTION 



Much recent work has raised the interest in fundamental 
questions on energy transport and heat flow at the nanoscopic 
scale where quantum effects become dominant.''^ Meso- 
scopic conductors are particularly suitable for investigations 
of nonequilibrium phenomena in quantum systems since elec- 
tron transport can be created, manipulated and detected thanks 
to the use of electric biases or thermal gradients applied to 
electronic terminals coupled to the mesoscopic sample. In this 
way, thermovoltages generated in response to a temperature 
difference have been observed in quantum point contacts,— 
quantum dotsii^ and ballistic microjunctions.^ These exper- 
iments can be successfully explained using linear-response 
theoretical models.^ Much more scarce are contributions that 
explore the nonlinear regime of transport. We have recently 
proposed in Ref. a general scattering theory of transport 
for mesoscopic thermoelectric effects in the weakly nonlin- 
ear regime taking into account screening interactions out of 
equilibrium. A natural development of our theory would aim 
at considering nonlinear effects in the heat current flowing 
through a mesoscopic system. This is the goal of the present 
work. 

Heat currents can be driven by electric voltages or tem- 
perature differences. In the former case, a heat transfer ac- 
companies an electric current, which in the limit of low volt- 
age gives rise to the Peltier effect; in the latter case, a ther- 
mal gradient produces a heat flux, which in the limit of small 
temperature differences is well described by the Fourier law. 
Nonlinear effects in mesoscopic systems demand the applica- 
tion of large driving forces across small distances, a require- 
ment within the scope of today's techniques." Nonlinearities 
have been predicted to cause thermal rectification effects 
and low temperature cooling. A deep understanding of non- 
linear effects is also needed for a careful assessment of the 
device performance of heat engines'^, heat pumps'^ '^ and 
multiterminal heat-to-electric cuiTent converters. Our the- 
ory is based on a voltage and temperature expansion around 
the equilibrium state. To lowest order, the linear response 
coefficients depend on the electrostatic potential at equilib- 
rium only. Higher order terms in the expansion define weakly 
nonlinear transport coefficients that are functions of the po- 
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FIG. 1. (Color online) Sketch of a dot system driven out of equi- 
librium by voltage biases Vi and V2 and thermal differences Ti and 
T2. Ti and r2 denote the reservoir-dot tunneling rates to the external 
reservoirs and Ed represents to the dot energy level. U is the internal 
dot potential illustrated as a shift of the dot bottom band. 



tential landscape out of equilibriumi^ This dependence can- 
not be neglected and is crucial to formulate a gauge-invariant 
theory.^^ We determine the screening potential up to the first 
order in the driving forces. In the isothermal case, charge 
pile-up processes are governed by particle injectivities. Im- 
portantly, these injectivities break the Onsager symmetry un- 
der reversal of an external magnetic field, which has been 
confirmed experimentally. ^''"^'^ In the isoelectric case, charge 
can be injected to or from the system depending on whether 
the caiiier's energy lies above or below the chemical poten- 
tial. Thus, we define entropic injectivities'^'^ that specify the 
charge pile-up in response to a pure thermal gradient. The 
electrostatic potential is then completely determined once the 
bare and the screening charges are calculated. As a result, the 
transmission function becomes a function of energy, voltage 
and temperature shifts.'" Below, we discuss the consequences 
of this important result for the nonlinear regime of heat trans- 
port. 

We examine two applications of our theory. First, we con- 
sider nonlinear contributions to the Peltier coefficient. Strik- 
ingly enough, there are very few works devoted to the Peltier 
effect beyond linear response. Exceptions are Ref. on 
metallic constrictions, Ref. 30 on quantum point contacts and 
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Ref.m] on bulk semiconductors. We consider a quantum dot 
attached to two terminals and calculate self-consistently the 
Peltier coefficient for higher-order currents. Interestingly, we 
find that the weakly nonlinear contributions are expressed in 
terms of a ratio difference that quantifies the relative impor- 
tance of the nonlinear conductances with respect to the linear 
ones. Additionally, we test our analytical results with numer- 
ical calculations of the full model. 

Second, we examine departures of the Wiedemann-Franz 
law out of equilibrium. We recall that this rule establishes 
a proportionality relationship between the linear thermal and 
electric conductances and that the proportionality factor, the 
Lorenz number, depends on universal physical constants only. 
Naturally, the Wiedemann-Franz law is not valid outside 
the Sommerfeld theory of metals (noninteracting electrons 
and linear-response regime), although it has been empirically 
found to be satisfied in bulk metals within a wide range of 
temperatures. Violations of this law in mesoscopic systems 
have been investigated in artificial Kondo impurities,^22iS in 
the fluctuating dynamics of open quantum dots,"*"* in single- 
electron transistors,^^ in strongly interacting dots coupled to 
ferromagnetic leads^'' and in double quantum dots.^^ Here, 
we discuss deviations from the Wiedemann-Franz law that 
emerge in the nonlinear regime of transport only. We find that 
these departures become maximal when the distance between 
the Fermi energy and the dot level is of the order of the level 
broadening because in that case the dot transmission energy 
dependence becomes strongest. 



II. THEORETICAL MODEL 

Our system consists of a mesoscopic conductor attached to 
multiple terminals a, /3 . . .. Each terminal is electrically and 
thermally biased as — Ep + eVa and 9 a = Tq, — T, respec- 
tively, iia being the electrochemical potential, Ep the Fermi 
energy, Ta the temperature of lead a and T the background 
temperature. Charge and heat transport is completely char- 
acterized by the scattering matrix Sap — SapiE, ell). Gener- 
ally, the scattering matrix depends on the carrier energy E and 
the electrostatic potential inside the sample [/.^'-^^ The inter- 
nal potential U = U{f, {Vy}, {0^}) is a function of the po- 
sition r and the set of applied voltages {Vy} and temperature 
gradients {(?-y}.'" Therefore, a complete calculation of Sap as 
a function of the full potential landscape U in the presence 
of interactions out of equilibrium seems an insurmountable 
task. We will then focus on the weakly nonlinear transport 
regime for which the transport coefficients and the system's 
response can be expressed in terms of quantities evaluated at 
equilibrium,^' 

We denote with Iq, and J7q the charge and heat currents, 
respectively, that flow from the leads toward the sample: 

/a = / dEA^p{E,eU)Jp{E), (1) 

Jc. = ^Y.J '^^(^ - ^^a)Ao.piE, eU)fp{E) , (2) 



where A^p = Tr[5ai3 - s^pSap] and fp{E) = 1/(1 + 
exp [{E — Ep — eVi3)/kBTp]) is the Fermi distribution func- 
tion in reservoir /3. A treatment of Eq. ^ in terms of interact- 
ing Green functions can be achieved following Ref. 38. Our 
discussion here is entirely based on the scattering approach. 
The sum over all heat flows is ^ • J7q = 

Eo^Jf - Eo^^cV^^ where 
{2/h)Y,f}! dEEAap{E,eU)fp{E) is the energy cur- 
rent. For time-independent driving forces, the energy 
current is conserved"*^ since unitarity of the scattering matrix 
{Y,a = 0) dictates that = 0. Heat fluxes thus satisfy 
the sum rule Y^ai'^a + ^Tq Vq) = 0. 

Equations ([T]l and (|2]l are exact within the scattering ap- 
proach. Now, in the weakly nonlinear transport regime charge 
and heat currents can be expanded around the equilibrium 
state (defined with Va —0 and 9^ — for all a) up to second 
order in powers of the driving fields Va and 0^ ■ 

la — ^ GapVp + ^ LapOp + ^ Gap-fVpVy 

P P Pi 

+ ^Lapy6p6-y + 2^Map^Vp6-y, (3) 

Pi Pi 

Ja = ^ RapVp + ^ KapOp + ^ Rap^VpVy 

P P Pi 

+ ^Kapy0p6^ + 2^Hap^Vpey. (4) 

Pi Pi 

The linear-response electric conductance is given by, 

Gap = ^J dEAap{E) [-dEf{E)] . (5) 

At very low temperature, the Sommerfeld expansion to lead- 
ing order in ksT/Ep yields the simple expression G ~ 
(2e^/ft.) Aap{EF)- The linear-response thermoelectric con- 
ductance is 

L-'P^J^j d^E^E ~ EF)Aap{E) [-dEf{E)] , (6) 

which reduces to L^p = 

{2e/hT){TT^k%/3)dEAap{E)\E=EF after the Sommer- 
feld expansion is applied. The linear-response heat current in 
Eq. (IDi is given by the electrothermal conductance 

Rc.p^^ J dE{E - EF)Aap [E] l-dEfiE)] , (7) 

and the thermal conductance 

Kap dE^^^^^AapiE) l-dEfiE)] , (8) 

where analytical expressions can be obtained for these 
conductances considering again the Sommerfeld expansion: 
Rap = i2e/h)in^klTy3)dEAapiE)\E=E^ and K^p = 
{2/h){7r'^k%T/3)AapiEF)- Note that the nondiagonal trans- 
port coefficients produce electric current from a temperature 
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difference (the thermoelectric conductance Lap) or heat cur- 
rent from a voltage difference (the electrothermal conduc- 
tance Rap)- Due to reciprocity, both effects are connected: 

Rap — TLap- 

The linear response coefficients in Eqs. (|5), dD, d?) and 
dHJ are evaluated at equilibrium and consequently Gap, Lap, 
Rap, and Kap are independent of the screening potential U . 

I 



On the contrary, the leading-order nonlinearities (see below) 
do depend on the applied electrical and thermal shifts through 
the electrostatic potential U . Indeed, this dependence is re- 
sponsible for rectification effects in both, electrical and heat 
currents. The nonlinear transport coefficients for transport of 
charge were obtained in Ref. [T^ We here find the nonlinear 
coefficients for the heat current: 



R, 



J dEdEf{E)Ya^Aap + SapAap - {E ^ Ep) 

(E-Ef)^ j fdAap , dAa 



fdA. 



KaPy = ^ / dEdEf{E) ^ 



H, 



80. 



803 



ap.^-J<iEdEnE)iE-EJ 



ap 



dAa 



\ deV-y deVp 
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■ Sp. 



7 ^'^P 

E — Ep dAap 



T deV. 
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T 

AaP 



op^ 



{E - Ef) 



dAap 
dE 



{E~EF)dAap , Aa^ 



T 



dE 



A 



aP 



dE T 

-{E~EF)dAap , AaP 



T 



(9a) 
(9b) 
(9c) 



Equations ( |9al l. ( |9b] l. ( |9c] l. are formally the main results of 
this work. Notably, the nonlinear responses depend not only 
on AaP but also on its change with variations of the set of 
electrical and thermal shifts {¥^,0^} due to the screening re- 
sponse of the system. 

A reasonable model for the interactions as response to the 
electrical and thermal biases considers only small deviations 
away from equilibrium. Hence, the internal potential is ex- 
panded as 



U ^U^^ + ^UaVa+^. 



(10) 



where Ua — {dU/dVa)eq and Za = {dU/d9a)eq are the char- 
acteristic potentials. These potential susceptibilities are mea- 
sures of the internal reaction of the conductor in response to 
an electrical and thermal shift applied to contact a. For defi- 
niteness, we next consider the case of a homogeneous poten- 
tial profile independent of the position, although we empha- 
size that the extension to inhomogeneous fields in our general 
model is straightforward. 

We first evaluate the total charge q of the conductor, which 
has two contributions: (;) (?barc corresponding to the bare 
charge injected from lead a and ( //) the screening charge de- 
noted with Qsci- The latter corresponds to the charge that 
builds up inside the conductor in response to the injected 
charges.^' The bare charge qhmc injected from lead a origi- 
nates from a voltage and a temperature imbalance in that ter- 
minal. Therefore, gbare is decomposed into the particle injec- 
tivity I'^iE) contribution^' and the entropic injectivity i^^{E) 



term:- 
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P 



dspa 



dE 



E ~ Ef t dspa 



T 



"Pa 



dE 



(11) 
(12) 



Importantly, the contribution to ^barc due to a temperature 
shift can be either positive or negative depending on whether 
the energy of the carriers is above or below Ep- This noto- 
rious feature is related with the fact that a heat addition or 
removal depends on whether the carrier energy E is larger or 
smaller than Epi^ Then, the factor {E — Ef)/T represents 
the entropy transfer associated to the additional carrier into the 
conductor We also note that in Eq. ( fT2] i we write Ef instead 
of the chemical potential /i(T), which is in general temper- 
ature dependent if the reservoir's charge density is assumed 
to be fixed.'*" However, in experimentally relevant situations 
one externally fixes the electrochemical potential using, e.g., 
voltage sources. Obviously, fj, ~ Ef in the limit of very low 
temperatures. Finally, the total accumulation or depletion bare 
charge imbalance becomes 



9bar( 



eJ2iKeVa + D^a0a), 



(13) 



where DP = - J dEiyP{E)dEf, and 

— J dEiya{E)dEf represent the integrated particle and 

entropic injectivities around the Fermi energy. 

The screening charge is calculated from the response of the 
internal potential, AC/ = U — U^q, away from the equilib- 
rium state C/eq. We consider the random phase approxima- 
tion, in which case the screening charge is proportional to 
the Lindhard function 11, gscr = e^HAU, which in the long 
wavelength limit becomes H = J dED{E)dE.f,'^^ with D = 
D{Ef) the conductor density of states. The set of equations 
for the characteristic potentials is closed using the Poisson 
equation. Thus, we relate the out-of-equilibrium net charge 
Sq = q - geq to AU = U - U^q through V^AC/ = -47rAg. 
By employing Eq. ( flOb and the fact that Va and Oa shifts are 
independent, we find a set of equations for the electrical and 
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thermal characteristic potentials: 



(14) 
(15) 



It is computationally useful to aply the WKB approxima- 
tion to the nonlinear coefficients Eqs. ( |9a] l. ( |9b] i and ( |9c] i. No- 
tice that this approach has the same range of validity than the 
long wavelength limit taken above. Therefore, we can make 

I 



the replacement S/5U — > —ed/dE in Eqs. ( |9a] i. ( |9b] i and 
( |9cl ) and the voltage and temperature derivatives, dg_^Aai3 and 
9v^Aq/3, are calculated once the characteristic potentials are 
known since 

de^Aap = ZjSAap/SU -ez^dsAap , (16) 
dv^Aap = Uy6Aap/6U -euydEAap . 



Thus, Eqs. 



and (I9ct become 



(9£; 



OA, 



9£: 

-{E~Ef) dA^p A^p 



iE-E,f4^^A^p 



dE 





9£; 
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T 9£; 



iE^EF)dAa,p , A 



(17a) 
(17b) 



ap 



T 



(17c) 



In the low temperature limit and for a two terminal device 
(a = 1, 2), the number of nonlinear conductances in Eq. (fTTj 
can be greatly reduced because the electrothermal and thermal 
conductances for contact a — 1 can be expressed in terms of 
the characteristic potentials ui 2, and zi 2 and the total trans- 
mission T{E) as 



^2 r 



^111 f 

Kni' 
Etiii '' 



h 3 



2zi 



dT{E) T{E) 
dE 



eT 



E — Ep 



e2 TT^klT^ 



zi- 



d^TjE) 1 OTjE) 
'dE^^''^^~dE~ 



E — Ep' 

(18a) 
(18b) 



E—Ep 

(18c) 



III. QUANTUM DOT 

In order to illustrate the general formalism described in 
the previous section, we now investigate the nonlinear heat 
transport for a paradigmatic mesoscopic system: an interact- 
ing quantum dot. Recent experimental findings support inter- 
esting nonlinear thermoelectric effects in quantum dots.'* We 
theoretically model the quantum dot system with a single lo- 
calized level with energy Ed- The dot is attached to two reser- 
voirs via tunneling barriers as shown in Fig. [T] We model 
such rates with energy independent constants Fi and F2, re- 
sulting in a total dot level broadening given by F = Fi + F2. 



We treat the electrostatic interaction within a mean-field de- 
scription. Under these considerations, the heat current from 
contact 1 is calculated from 

Ji^lJ dE{E-EF-eVi)r{E) [h{E) - ^(-B)] , (19) 



where 



T 



4Fir2 



{E-Ed~ eUy + F2 



(20) 



is the corresponding Breit-Wigner transmission Une shape for 
the dot level as a function of the internal potential U . The dot 
chai-ge reads, 



ttJ {E 



T^h{E)+T2f2{E) 



Ed - eUf + F2 



(21) 



Within a mean-field treatment of interactions, we solve the 
discrete version of the Poisson equation by introducing a ge- 
ometrical capacitance C which connects electrically the dot 
to an external gate terminal Vg controlling the level position. 
Then, U is determined from 



5qd = qd- qT = C{U - Vg) 



(22) 



where 5qd represents the charge excess due to electrical and 
thermal applied biases while cf^ is the equilibrium charge cal- 
culated from Eq. ( 1211 1 by setting /i — f2 — f ■ 

Once the internal dot potential U is obtained, Eq.(fT9]l can 
be integrated and yields 



.^2 



(23) 
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where ^' is the digamma function with arguments ^1(2) = 

/3i(2)[^d ~ iT ~ /xi(2)] where ^1(2) = l/fesT^z). Im- 
portantly, the dot level is renormalized by the interactions: 
Ea = Ed + eU. 

Equation |23] is the exact expression for the heat current 
when Coulomb interaction is considered within a mean-field 
description. It is to be compared to the voltage-temperature 
expansion obtained from the nonlinear heat transport formal- 
ism discussed in Sec. II. To do so, we calculate the heat trans- 
port conductances in terms of the characteristic potentials and 
expand Eq. (ISTT l to leading order in Va, Oa and U. We find 



Sq^ = e^DfVi + e^ZJf F2 + eD'^Oi 



e^DU, (24) 



where 



DP 



Df. 



n 



dE 



dE 



dEf{E) 



{E - EdY ■ 
E — Ep 



_ p2 ' 

dEf{E) 



T [E- EdY + r2 



(25) 



(26) 



Notice that Eqs. ( |25l l and ( |26] l are the integrated particle and 
entropic injectivities within a Breit-Wigner representation of 
the dot scattering matrix. Using Eq. (l22t and (|24] | we find the 
dot internal potential 



U 



e^DlVi + e^DlV2 + eDlOi + eDf 6I2 + CVg 



C + e^D 

from which the characteristic potentials follow. 



(27) 



Ul(2) 



1(2) 



c 



eD 



1(2) 



C + e^D' 



C + e^D 



C + e^D 



Zl(2) 

(28) 

g = 1 because of gauge 



Note that while ui + U2 + u 

invariance^ such a sum rule is not satisfied by the z poten 
tials . 
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FIG. 2. (Color online) Heat cuiTent J (in units of Jo = 2VxV2lh) 
for a quantum dot system with /ii = Ef -I- eV , fj.2 = Ef, 
61 — T + 9, and 62 — T for three different values of the tunneling 
asymmetry rj and Ed = lOF with T = O.SF. (a) J-V characteris- 
tics (full lines) for 6* = along with the leading-order nonlinearity 
J' ~ RiiV + RiiiV^ (dashed lines) for the 77 values indicated in the 
right panel. Upper inset: we show with solid lines the self-consistent 
potential [/ as a function of eV/T. Dotted lines correspond to the 
approximated values U ~ uV = Fi/F to leading order in an ex- 
pansion in powers of V. Lower inset: temperature dependence of 
(b) J'-9 characteristics (full lines) for V — along with the 
leading-order nonlinearity J ~ KnO + KmO^ (dashed lines) for 
the indicated values of 77. Upper inset: we show with solid lines the 
self-consistent potential f/ as a function of the thermal shift. Lower 
inset: temperature dependence of the nonlinear thermal conductance 



IV. HEAT RECTIFICATION EFFECTS 

Rectification effects occur in systems where the functional 
dependence of current versus the driving field (voltage or tem- 
perature shift) departs from being linear due to quadratic (non- 
linear) transport responses. We now investigate nonlinearities 
in the heat flux-voltage and heat-temperatures characteris- 
tics of a quantum dot. We consider the charge neutral limit 
(C = 0), in which case strong interactions renormalize the 
energy level to maintain a fixed charge inside the dot. This 
limit is relevant in many experimental situations. 

In the isothermal case (Ti = T2), we take V\ — V and 
V2 = 0. Then, to leading order in voltage the character- 
istic potential becomes u = dU/dV = (1 -I- 77)/2 with 
r] = (Fi — r2)/r the tunneling asymmetry.-- In the upper 
inset of Fig. |2a) we show the exact dot potential U obtained 
from a numerical, self-consistent calculation of Eq. (l2ll com- 
pared to its leading-order value U ~ uV . We observe that 
the approximation is rather good even for voltages larger than 
the level broadening F. The heat current ^7 = J7i is shown in 



the main part of Fig. |2]as a function of voltage. At low volt- 
ages, the relation is linear but is not clearly visible since the 
quadratic-in-y term largely dominates (dashed lines). We re- 
call that at low temperatures J ~ RmV'^ represents a Joule 
heating term, which is the main source of heat release in a 
voltage-biased terminal. In the lower inset of Fig. |2la) we 
plot the nonlinear electrothermal conductance as a function 
of the background temperature. Interestingly, i?iii increases 
in magnitude more quickly as the tunneling asymmetry in- 
creases. 

In the isoelectric case (/ii — /12), we take Oi — and 
02 = 0. The characteristic potential determined from the en- 
tropic injectivities is z ^ dU/de ^ D1/[2e{DP+D^)]. Thus, 
the dot potential reduces to U — zO to leading order in 0. We 
compare in the upper inset of Fig. |2lb) the validity of this 
approximation with the exact calculation of U . We find that 
for thermal shifts much smaller than F, the relation between 
U and 9 stays linear. In the main part of Fig. |2lb) we de- 
pict the heat current as function of for various values of the 
tunneling asymmetry. Since we select small values of 0, the 
agreement between the weakly nonlinear expressions and the 
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full calculation is excellent. In the lower inset of Fig. |2b) we 
show the nonlinear thermal conductance Km as a function 
of the background temperature T. Similarly to the isothermal 
case, the rectifying thermal conductance increases with higher 
values of T. 



V. NONLINEAR PELTIER EFFECT 

The Peltier effect accounts for the generated heat in a given 
terminal caused by the flow of a unit charge through the con- 
ductor, keeping all terminals at the same temperature. More 
precisely, the two-terminal Peltier coefficient is defined at lin- 
ear response as the ratio between heat and electrical currents. 



n„.f 



Rii 



(29) 



when 9y = Q for all 7. Thus, the Peltier coefficient can be 
viewed as the analogue to the thermopower Sq (a voltage gen- 
eration in response to a thermal gradient under the condition 
of zero net current). In fact, both coefficients are connected 
via the Kelvin relation XIq — TSq. The Peltier effect informs 
us about how efficiently electrical currents are converted into 
heat currents whereas the Seebeck effect quantifies the trans- 
formation of waste heat into useful electricity. 

Our goal here is to generalize Eq. (|29^ to the nonlinear case. 
We consider the following expression. 



n = — 



(30) 



To calculate the Peltier coefficient in the weakly nonlinear 
regime we employ the second order expansion in terms of 
electrical and temperature gradients for both the electrical and 
heat currents. For a two-terminal device an isothermal electri- 
cal current 2i = I is driven by a bias voltage drop V, 



I = GiiV + GmV^ + OiV"^) . 



(31) 



Inverting this expression, we substitute the voltage V = 
X/Gii - I^Giii/Gii + 0{I^) H into the heat cuiTent, 



(32) 



yielding 



Ji 



Rii 



1 

Gil 



^111 Gil 



Gil Gil V -^11 ^11 
We insert Eq. (l33T l in Eq. (|30] | and find 

Gill 



n = nn 



1 



1 (Rill 

Gil 



Rii Gil 



+ 0(13) ^ (33) 



(34) 



where we have set 11 = Hi. Now, the second term in the 
r.h.s of Eq. (|34] i measures the deviations of 11 from the linear- 
response value, Hq. Interestingly enough, nonlinear contribu- 
tions are dictated by a ratio between the nonlinear to linear 
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FIG. 3. (Color online) (a) Linear- response Peltier coefficient Ho for 
a current biased quantum dot at Ed ~ lOT. We show both the low 
and high temperature limits, (b) Full Peltier coefficient 11 including 
contributions beyond linear response for three different background 
temperatures foi rj = 0.5. We show with colored dotted lines the 
leading-order expansion calculated from Eq. i34\ . 



thermal and electrical conductances. In other words, the con- 
version efficiency for electric currents into heat flow is given 
by the relative strength of the nonlinear conductances to the 
linear ones and by the difference between the heat and the 
electrical properties of the conductor 

We show in Fig.[3la) the linear Peltier coefficient for a two- 
terminal quantum dot as described in Sec. II taking Gm — 
(eVM (1 - 2wi)aBr(-E) Ib^b^ As the background T in- 
creases, XIq quickly enhances following a law and then sat- 
urates for temperatures much larger than the resonance broad- 
ening r. The first fact can be understood from a Sommerfeld 
expansion of i?ii and Gn. The former goes as whereas 
the latter is independent of temperature to leading order in T. 
Finally, the saturation value can be deduced by replacing the 
Lorentzian shape of the Breit-Wigner resonance with a delta- 
like peak, an approximation valid in the high temperature 
regime (T > T). Then, one finds Hq = -{Eo - E - F)/e. 
Note that this expression obeys the Kelvin relation XIq — TSq 
since the linear-response thermopower for the same system is 
^ -{Eo- E ~ F)/eT at high temperatures. 1° 

In the nonlinear regime, 11 depends on the nonlinear re- 
sponses according to Eq. ( [34l i and a self-consistent calculation 
is thus needed. We depict in Fig.[3jb) the exact value of 11 ob- 
tained from a numerical calculation of Eqs. ([T), and (|2) in the 
mean-field approach for interactions (full lines), together with 
the weakly nonlinear result expressed by Eq. ( |34t . To leading 
order in I/Iq, where Iq = {2e'^ /h)riT2/r, the Peltier coef- 
ficient shows a linear dependence of the current, in agreement 
with Eq. ( l34l i. Furthermore, 11 increases for increasing I but 
only in the low current regime. At higher currents, highly non- 
linear terms start to contribute and deviations from the linear 
regime are observed. 
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FIG. 4. (Color online) Nonlinear departures of the Wiedemann- 
Franz law as expressed by Eq. (139) . Nonlinear transport responses 
are calculated for a two-terminal quantum dot with voltage bias 
V = 0.25r, background temperature T = O.IF and different values 
of thermal shift as a function of the dot level position Ed- 



VI. WIEDEMANN-FRANZ LAW 



. Clearly, Eq. (I35l l reduces to Eq. (I38l l if nonlinear terms are 
neglected. Departures from the Wiedemann-Franz rule occur- 
ring in the nonlinear regime of transport can thus be quantified 
with the following magnitude: 



A-Ao ^111^ Gill, 



Gil 



(39) 



Interestingly, deviations from the Lorenz number are given 
by the difference between the ratios of nonlinear to linear ther- 
mal and electrical conductances. 

We calculate Eq. (3% for the quantum dot model described 
in Sec. Ill for V = 0.25r and T = O.Oir and different val- 
ues of 6. We show our results in Fig. |4]as a function of the 
dot level position Ed- We find that departures are stronger 
for increasing temperature differences since for increasing 9 
the heat transport becomes more nonlinear. Moreover, when 
Ed Hes around F above or below the Fermi energy departures 
are more apparent because the transmission function changes 
rapidly for energies around those points. 



The Wiedemann-Franz law for metals establishes that the 
ratio between the linear-response thermal (K) and electri- 
cal (G) conductances is proportional to the temperature (T) 
with a proportionality constant given by the Lorenz number, 
Aq = (7r^/3)(fcs/e)^. This relation holds as long as the Som- 
merfeld expansion is valid. More concisely, the Wiedemann- 
Franz law for a mesoscopic system states that 



K 
TG 



An 



(35) 



which is obviously satisfied in the linear regime since K ~ 
{2/3h)T:^klTT{EF) and G ~ {2e^ /h)r{EF)- 

To generalize this result to the nonlinear regime and to test 
whether the Wiedemann-Franz law is still valid or departs 
from the Lorenz value we consider the Wiedemann-Franz ra- 
tio 



Ar 



{V=0} 



T{IJV) 



(36) 



{9=0} 



For illustrative purposes we have considered again the case 
of a conductor attached to two terminals, where only one of 
the contacts is heated or cooled with thermal difference 6 and 
electrically biased with voltage V; e.g., di ~ 9, 62 ~ 0, Vi ~ 
V, and V2 = 0. Then, 



A = 



1 Kn+Kni9- 
TGn+GinV 



(37) 



Keeping only terms linear in 9 and V, the nonlinear analogue 
to the Wiedemann-Franz ratio is 



A^l 
T 



Gil 



KnGinV , Knit 



ri2 



Gil 



(38) 



VII. CONCLUSIONS 

We have investigated the heat flux in a mesoscopic system 
beyond linear response. Our scattering approach includes in- 
teractions to take into account the system's screening response 
to voltage and temperature shifts in the reservoirs. We have 
found exact expressions for the leading-order nonlinearities 
in a voltage-temperature expansion of the heat current and 
have applied our general theory to the study of the heat rec- 
tification properties of an interacting quantum dot attached to 
two terminals. We have found that the Peltier coefficient ac- 
quires higher-order contributions at moderate electric currents 
and that these contributions are stronger for increasing back- 
ground temperature. Furthermore, we have investigated de- 
partures from the Wiedemann-Franz law in the weakly nonlin- 
ear regime of transport and have found that the strongest devi- 
ations occur when charge fluctuations dominate, i.e., when the 
distance between the dot level and the leads' Fermi energy is 
of the order of the resonance broadening. We believe that our 
results are of fundamental importance for the understanding 
of basic nonUnear heating effects in quantum conductors. 
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